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Abstract 

We propose a model that provides a unified description of nuclear equation 
of state and fragmentations. The equation of state is evaluated in Bragg- 
Williams as well as in Bethe-Peierls approximations and compared with that 
in the mean field theory with Skyrme interactions. The model shows a liquid- 
gas type phase transition. The nuclear fragment distributions are studied for 
different densities at finite temperatures. Power law behavior for fragments 
is observed at critical point. The study of fragment distribution and the 
second moment 52 shows that the thermal critical point coincides with the 
percolation point at the critical density. High temperature behavior of the 
model shows characteristics of chemical equilibrium. 
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I. INTRODUCTION 



The phenomenon of nuclear matter fragmenting into various pieces can be studied in 
heavy ion collisions. This has been an area of much activity, both in experiments and in 
theory. Curtin, Toki and Scott 0] pointed out that at some incident energies excited matter 
that is formed in heavy ion collisions will pass through a liquid-gas phase transition stage 
and if fragments are formed at this stage, it may show characteristics of this transition. A 
study of nuclear matter with a Skyrme interaction was made in || where it was shown that 
in mean field theory there is a phase transition as in Van der Waals gas. In nuclear physics, 
phase-transition, if it indeed happens, is a transient phenomenon and it is not clear what the 
mean field theory predicts for observables that can actually be measured. In experiments 
the most readily accessible observable is the cross section of nuclear fragments, or yield Y(A) 
vs. A where A is the mass number or charge of a fragment. Of course mean field theories can 
not make a prediction about fragment distribution thus it falls short of directly providing 
results with which one could confront the data. 

Bauer 0] and Campi j| used the percolation model || to calculate fragment distribution 
H|7]]. There are two varieties of percolation models: the bond and the site percolation. In a 
bond percolation model each site of the lattice is occupied by a nucleon. That is, the number 
of nucleons equals the number of lattice sites N. The bonds between nearest neighbors are 
broken with a probability of 1 — p where p is the probability that the bond is unbroken. 
Nucleons which are connected through unbroken bonds form a cluster. In a site percolation 
model, each site is occupied randomly with probability p < 1. The number of nucleons, n, 
is usually less or equal to the number of the lattice sites N. Nucleons occupying nearest 
neighbor sites are considered to be the part of one cluster. In both bond and site percolation 
models there is a critical value p c of p at which an infinite cluster starts to emerge. For very 
large lattice p c is independent of N. The probability that a given site belongs to this 
infinite cluster is zero for p < p c , and grows from to 1 for p > p c . In general, the cluster 
distributions in percolation models are very similar to the mass distribution observed in 
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heavy ion collisions 

The percolation model is quite different from mean field theories that we first alluded 
to. There is just one parameter in a percolation model. In a bond percolation model the 
number of sites equals the number of nucleons. One might regard that the probability p is 
a function of temperature. In that case there is no parameter corresponding to the volume 
or the pressure of the system. If instead one takes the site percolation model, the number 
of lattice sites N is usually larger than n. One may now associate the parameter p with the 
volume. There is, then, no reference to temperature or pressure. Thus the simple percolation 
model can not describe the thermodynamic aspect of nuclear fragmentation. 

Our present work started with the desire to have a model that can both describe the 
equation of state and the fragmentation of finite nuclear systems. The quintessence of this 
model is the lattice gas model ||. The model leads to a (P,V,T) diagram, and it is also 
linked with the percolation model in an obvious way. An interesting feature is that it not 
only leads to a liquid-gas type of phase transition, but also encompasses the percolation 
transition. 

Much work has been published on the subject of liquid-gas phase transition in nuclei. See 
Refs. [P|-[TT[| for description of some early works and [|T2J for review. There are microscopic 



models which are proposed to calculate multifragmentation [|l3| -|l7ll- Other phenomenological 
models can be found in ||18|| . The pioneering work by the Purdue group |1| was a big impetus 
for this subject. 



II. THE MODEL 

We consider the participant zone in collisions of two heavy ions and conjecture that 
because of nucleon-nucleon collisions the system reaches thermal equilibrium. We assume 
that the system is adequately described by classical statistical mechanics. In that case the 
canonical partition function of an n-particle system can be written in a separable form: 
Z can oc Z p (can)Z r (can) where Z p (can) is given by 
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Z p (can) = / exp[-/3^ Pl 2 /2m]rf 3 Pl d 3 p n (2.1) 

J 1 

where (3 is the inverse temperature. The other part of the partition function is 

Z r (can) oc / exp[— /^T^ v(r ij )]d i ri....d i r n (2.2) 

i<j 

Here t> (r^) is the potential between particles i and j. We approximate the configuration part 
of the partition function by the partition function of the lattice gas model. In this model 
(as in a site percolation model) each lattice site can be occupied by at most one nucleon. 
The number of lattice sites, N, gives a measure of the volume of the system and is usually 
larger than n. When the nucleus is in the ground state we have n = N. Thus our model 
is limited to normal nuclear volume or higher. Becuase cluster formation presumably takes 
place at a volume much larger than normal volume, this restriction may not be a debilitating 
factor. In contrast with the percolation model, the lattice gas model includes interactions. 
If two nearest neighboring sites are filled by nucleons they will interact and the energy of 
interaction is denoted by — e. The nearest neighbor interaction simulates short range nature 
of nucleon-nucleon interaction. 

We can now write down the canonical partition function for our model for a fixed n and 
N. The Z p (can) is simple and does not depend on the volume, and we just need to calculate 
the Z r (can). Let N nn be the number of nn bonds in a specific lattice configuration, the 
energy carried by these bonds is then —eN nn . Then Z r (can) is given by 

Z r (can) = ]T g(N, n, N nn )e^ N ™ (2.3) 

Nnn 

where g(N,n, N nn ) is a degeneracy factor satisfying 

ATI 

P^ NM = W^- {2A) 

An exact evaluation of Z r (can) given in (|2.3| ) is difficult. Thus we will use approximate 
means. 
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III. EQUATION OF STATE 



A. Bragg- Williams approximation 

The Bragg- Williams approximation is an easy and quick calculation but is not expected 
to be accurate. The results are, however, transparent and analytically demonstrate Van der 
Waals type behavior. In this approximation the number of nn bonds N nn is taken to be 
given and fixed, when N and n are specified. If one site is definitely occupied, the number of 
its 7 neighbors that are occupied is on the average jn/N. In our three dimensional simple 
cube lattice 7 is 6 except at the boundary. Since there are n nucleons in the system, the 
number of nn bonds is then r yn 2 /2N where we have ensured that each bond is only counted 
once and assumed that both n and N are large so that boundary effects can be neglected. 
In Bragg- Williams approximation the canonical partition function is then 

z ^ = w^ e¥,li - (31) 

The equation of state can be calculated by utilizing P = kT(d\n Z(can) / OV)t = 
kT(d\n Z r (can) /dV) T since Z p (can) does not have any V dependence. Here P is the pres- 
sure and V is the volume given by V = a 3 N. A representative value of a 3 would be 
a 3 = l/po = 6.25 fm 3 where po is the normal nuclear density. The normal nuclear volume is 
Vo = a 3 n. Using Stirling's approximation for N\, n\ and (N — n)\ one can show 

„ kT , N 1 . n ,„ ^ x 

^ = ^F^-^<^ 2 - < 3 ' 2 > 

Using n/N = V /V = pj p we finally get 

This equation of state has the same qualitative behavior as the Van der Waals gas. For one 
mole of gas, the Van der Waals equation of state is 

NfykT a 



The lattice gas pressure goes to infinity as V approaches Vq. The Van der Waals gas pressure 
goes to infinity as V is squeezed to the value b. For large V both the equations of state 
approach the perfect gas limit. The critical point can be readily determined analytically 
from (|3.3|). By setting dP/dp = d 2 P/dp 2 = at the critical point we obtain pc = 0.5po 
and kTc = 7^/4, respectively. It is also straightforward to show that in the Bragg- Williams 
approximation PcVc/RTc = 2 In 2 — 1 = 0.386 for the lattice gas. The corresponding 
number for the Van der Waals gas is 0.375. 

The Bragg- Williams approximation is considered to be crude and one may wonder if 
the lattice gas model would indeed lead to a liquid-gas type phase transition when a better 
mean-field calculation is done. In the next subsection we try an improved approximation. 



B. Bethe-Peierls approximation 

We now try to do a better mean field calculation and use what is called the Bethe- 
Peierls approximation in the Ising model. For the Ising model, the order parameter can be 
computed without having to calculate the partition function but for the equation of state that 
we wish to calculate we will need to obtain the partition function. Here we consider a grand 
canonical ensemble. To explain the methodology we refer to Fig. 1 where, for simplicity, a 
two dimensional square lattice is shown. We break up the lattices into N/ (7 + 1) blocks, 
each of which contains 7 + 1 sites. The interactions within each block are treated exactly 
while the interactions between different blocks are taken into account approximately. The 
local correlations are taken into account in this approximation, and it is an improvement 
over Bragg- Williams approximation. The grand partition function can be written as the 
product of the grand partition functions of the N/i^y + 1) blocks: 

N 

Z gr = z gr (b\ock l),2 3r (block 2) • • • z 9r (block ). (3.5) 

7 + 1 

The grand partition function of the block labeled by 1, 2, 3, 7 and 5 can be written as 

7+1 

v = EE^^ EoW (3-6) 

k=0 a 
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Here A' plays the role of chemical potential. The energy E a (k) consists of two parts; 
the kinetic energy (Yli=i P 2 i / an d the potential energy for nearest neighbor interac- 
tion. Integration over the kinetic energy part can be immediately done to give a factor 
of [(2irm//3)^ 2 ] k = e qk . We now define A = A' + q and also divide the the right hand side of 
the above equation into two parts: 



z gr 



£ h\ e Xk +m + e X £ h\ e Xk+me+e)_ ( 3 7 ) 

k — \ / k — \ / 



The first part is proportional to the probability that site 5 is unoccupied and its neighboring 
sites 1 • • - 7 are fc-fold occupied where k goes from to 7. The second part is proportional 
to the probability that site 5 is occupied where an extra amount of energy —ke is included 
when k of the nearest neighbors are also occupied. The extra factor exp[/5fce] takes into 
account the interaction between different blocks. If we left this factor out the nucleon at 
site 1 (Fig. 1) would only interact with the one at site 5. 

Because of self-consistency condition the average occupation probability at every site 
must be the same. The average occupation at site 5 is 

^(5) = L = I (3-8) 



Zgr V 



where 



z gr = (1 + e x+pl V + e A (l + e x+pe+pl y (3.9) 

is obtained from ( |3.7|) by summing over k. The average number of particles in all the sites 
neighboring site 5 is 

7 y ( -y\ faXk+fikl 1 A y h\ u \k+/3ke+/3ke 

P(i) = =^ =^ • (3.10) 

i=l z gr 

The self-consistency condition then implies 

e A (1 _|_ e \+pe+Pe\"t _ e^+Pen _|_ e X+pe\-y-l _|_ e 2X+P{e+e)/^ + e \+f3(e+I)y~l _ (3.11) 

The two unknowns e A and e can be solved from (|3.8|) , ( |3.9|) and ( |3.11|) . 



Divide both sides of (|3TTT| ) by e A (l + e ^(^))^-\ we obtain 

1+e A +/3? V -l 



1 = e 



/3? 



1 _)_ e A+/3e+/3? 



(3.12) 



and (|3.8|) can be rewritten as 



Use (gTp , (|3~T3l) becomes 



iV (1 + Ae^) 7 

rT ~ + e A (l + e A +^+^)^' 



— - 1 

n ~ + e A+/3? 7 /( 7 -l) • 



(3.13) 



(3.14) 



We can now write e in terms of e 



n 



-#=7/(7-1) 



Combining 



iV - n 

(|3.9| ) and (|3.15 ) we obtain 



(3.15) 



1 

x = - 
2 



N-2n ,N-2n„ 

1_ 4 / f ^2 _|_ 4(_ ^g/3e 



iV-n 



N — n 



'N - n' 



(3.16) 



where 



x 



(3.17) 



The values of e and e A can now be found from (|3.15|) - (|3.17|) . We note here that the results 
depend on the ratios of n/N and e/kT only. 

Let us now go back to the partition function for the lattice given in ( |3.5| ) where it is 
written as a product of the partition functions of the N/ (7 + 1) blocks. If we simply use the 
partition function for each block given in ( |3.9| ) we count twice the binding energies between 
neighboring sites in different blocks. For example, the binding energy between 1 and 6 (Fig. 
1) is included in z gr (block 1) through e, and it is included again in z gr (block 2). We note 
that on the average there are n/N particles at each site, and each block has 7 peripheral 
sites. Thus, when we evaluate the partition function for the lattice, the partition function 
for each block given in ( |3.9|) should be corrected with following factor: 



correction = e -^ n / N . (3.18) 

We can now use PV/kT = \nZ gr , V = N/p and \nZ gr = N/(j + 1) lnz gr to obtain 

P = p kT x — j-j- In z gr . (3.19) 

Here z gr is now understood as the product of a block partition function given in ( |3.9| ) and 
the correction factor given in ( |3.18| ). 

The equations of state calculated in Bragg- Williams approximation and Bethe-Peierls 
are compared in Fig. 2 in P — V diagrams at different temperature. It is seen that in high 
temperature limit both give the same results, and they begin to differ at low temperatures. 



C. Mean field theory with Skyrme interaction 

We take a Skyrme interaction with potential energy density given by 

V(P) = yPo(f ) 2 + -^TPo(f ) CT+1 (3.20) 
2 po a + 1 po 

where A = —356 MeV, B = 303 MeV and a— 7/6. This interaction produces a saturation 
density of 0.16 fm 3 and binding energy of 16 MeV per particle. 

The pressure consists of two parts: one originating from the interaction and the other 
from the kinetic energy. The pressure produced by the interaction is given by 



+ aB ( p y 
2 po a + Vpo } 



p. (3.21) 



The kinetic pressure is calculated numerically from a Fermi gas model at finite temperature. 

Similar to Fig. 2 one can draw P against V/Vq (= po/ p) at various T. A comparison 
of the equation of state in the mean field theory with Skyrme interaction and in the lattice 
gas is shown in Fig. 3. In both cases we see characteristics of liquid-gas phase transition. 
In the calculation we used e=9 MeV and po = 0.16 fm~ 3 . 
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IV. FRAGMENT DISTRIBUTIONS 



Nuclear fragmentations are described by the formation of clusters in our model. To gen- 
erate fragements we need to simulate the lattice configuration and momenta of particles, and 
to determine if neighboring particles belong to the same cluster. The lattice configuration is 
generated according to the ( |2.2j ), while the momenta are generated from Maxwell-Boltzmann 
factor given in ( |2.1| ). These two processes are independent of each other since the partition 
function can be factorized. To generate the configuration we start from an empty lattice and 
put the first nucleon at random. Once this has been put in the 7 boxes that are immediate 
neighbours are assigned a probability oc exp[/3e] whereas all other boxes have probability 
proportional to unity. The next nucleon is then put in according to this probability dis- 
tribution. If at an intermediate state there are m empty boxes we assign to each of these 
boxes a probability proportional to exp[g/5e] where q is the number of nearest neighbors 
that are already filled up. The next filling is then done according to this distribution. The 
difference from the computer simulation that would be done in a site percolation model is 
the Boltzmann factor. Having done the configuration space sampling we then assign each 
nucleon a momentum according to Maxwell-Boltzmann distribution. 

In a site percolation model two neighboring particles always form a cluster. The for- 
mation of a cluster in our model depends on the interactions of neighboring particles and 
their relative kinetic energy. We adopt the following physical criterion for determining the 
formation of a cluster. Two neighboring nucleons belong to the same cluster if the following 
condition is satisfied 

p^/2/x-e<0 (4.1) 

Here p r is the relative momentum between the two nucleons and [i is the reduced mass. 
Except at very low temperatures, the frequency with which two nucleons appear at neigh- 
boring sites depends mostly on density. The probability of p^/2/i exceeding the value e in- 
creases with temperature since the momenta of each nucleon is obtained from Monte-Carlo 
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sampling of Maxwell-Boltzmann distribution. Hence the probability that two nucleons are 
bonded decreases with increasing temperature and the system becomes less compact at 
higher temperature as it should. A different parametrisation used in |20] leads to similar 
effects. 

We note that when each particle obeys the Maxwell-Boltzmann distribution, the dis- 
tribution of relative momentum between two particles is also Maxwell-Boltzmann , i.e., 
P(p r ) = l/(27r/ifcT) 3//2 exp[— p^/2/i&T]. We can then write down a formula for bonding 
probability which is temperature-dependent: 



Coniglio and Klein [p2| used a different parametrisation for the bonding prbability. They 
used 

p= l-exp[-/3e/2]. (4.3) 

This was mathematically devised so that the thermal critical point would also be a percola- 
tion point, a feature, as we shall see, is also present in our parametrisation. A comparison 
of the two formulae is presented in Fig. 4. An example of fragment distribution obtained in 
our simulation is shown in Fig. 5. One could decipher from this figure that at T = 0.5Tc the 
system percolates; at T = l.bTc and T = 2.0Tc there is no percolation and that percolation 
sets in around T = Tc- 



V. THERMAL CRITICAL POINT AND PERCOLATION POINT 

Let p c be the value of p at the percolation point, i.e., the point when an infinite cluster 
just appears. The cluster size distribution for a general p can be parametrized as |J; 

Y(A,p) = A^f[(p- Pc )A°}. (5.1) 

The scaling function f(x) can be determined by computer simulations or experiments. In 
many cases, such as the Fisher model P5| , f(x) is an exponential function. At the percolation 
point, the fragment distribution obeys a simple power law 
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Y(A, Pc )ocA- T (5.2) 

At the percolation point the fluctuation is the maximum. This means that if in the neigh- 
borhood of the percolation point the yield Y(A,p) is fitted by a power law, the exponent r 
will be a minimum at the percolation point. Experimental data are often fitted to a power 



law and a minimum in r is searched to ascertain the percolation point |H[ . 

The second moment is defined as 5*2 = A 2 n(A)/n where n(A) is the number of clusters 
with A nucleons and the sum excludes the largest cluster (defined henceforth as A max ). In the 
thermodynamic limit S2 diverges at the percolation point and obeys a power law distribution 

S 2 oc \p-Pc\~ 1 - (5.3) 

As mentioned, the order parameter given by lim^oo A max /n is zero for p < p c in and 
increases with p for p > p c . Above the percolation point the order parameter is given by 

A m ax/n oc \p~Pcf- (5.4) 

There are two independent critical exponents in most statistical models and percolation 
models. One can readily show that following relations exist. 

3 — t , 

7 = ; (5-5) 

a 

P = — (5.6) 
a 

and 

r = 2 + -^-. (5.7) 

13 + 1 

For a finite system the second moment S% goes through a maximum (instead of infinity), 
and r goes through a minimum at the percolation point. Thus r and 5*2 can be used to 
identify percolation point in a finite system. 

In our model there is, first of all, a thermal critical point which is obtained at pc = 0.5po 



and Tq = 1.1275e [[24j]. This more exact value of Tq = 1.1275e is smaller than the value 
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of 1.5e that we obtained from Bragg- Williams approximation. But, in addition, there is 
a continuous range of percolation points. Provided the density is higher than a minimum 
value (ss .3po) a percolation point will be reached at a certain temperature. The higher the 
density, the higher the temperature at which percolation sets in. If we identify the point 
at which a maximum in S 2 or a minimum in r is achieved as the percolation point, then 
we see from figs. 6 and 7 that at higher density percolation point is reached at a higher 
temperature. One can now ask the question: if a minimum in r is seen in experiments as in 
ref [19], one could interpret that as an indication that the percolation point is reached but 
has it got any relevance with the thermal critical point ? As figure 7 shows if the freeze-out 
density is 0.5p then the percolation point is reached at the thermal critical point. It also 
follows that if the freeze-out density is close to 0.5p then a minimum in r is obtained when 
the temperature of the dissociating system is close to the critical temperature. 

We therefore have this remarkable result which was not a priori imposed. In our model 
the thermal critical point is also a percolation point. Numerically this result can be explained 
by noting that at the critical temperature the value of p in our model differs little from 
Coniglio-Klein parametrisation. At density p/p =l, the system begins to percolate at T = 
1.477c in our model and at T = 1.557c in Coniglio-Klein model. 

VI. MULTIPLICITY AS A VARIABLE 

The implicit thinking in much of what is presented above is that in nuclear collisions 
matter is compressed, heated up and we can talk of a freeze-out density and temperature at 
the time of dissociation. Both the temperature and the density are, however, not directly 
measurable. In the past values of |po an d lower have been used for the freeze-out density. It is 
not known accurately. Given that the most easily measurable quantities are the multiplicity 
(number of charged particles),^ and A max , we ask: can we determine the freeze-out density 
from the multiplicity dependence of £2 and A max ? In Figs. 8 and 9 we plot the S 2 and the 
A m ax as functions of the multiplicity at different densities. In the calculations, at a given 
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density we took a sufficient large range of temperatures such that we cover the full range of 
multiplicity. For a given multiplicity both S2 and A max can vary from one event to another. 
What is plotted is the average for a given multiplicity. From Figs. 8 and 9 we can easily see 
the changes of 5*2 and A max when the density is changed from 0.3po to 0.5po- The differences 
are however rather small between density 0.5po to l.Opo- Remembering that there will always 
be uncertainties in experimental data due to contamination from pre-equilibrium particles, 
spectators etc., we conclude that it is difficult from these observations alone to determine 
the freeze-out density accurately. Some other variables might better differentiate between 
different freeze-out densities. 



VII. BEAM ENERGY AS A VARIABLE 

Fig. 6 shows that it may be possible that freeze-out density can be determined from the 
temperature dependence of r. The temperature is not a direct observable in experiments 
although it has often been deduced indirectly from other data, notably from slopes of inclu- 
sive cross-sections. Here we will first try to deduce the temperature from a simple theory so 
that the temperature is given once the beam energy is given. We consider the experimental 
setup of [20|. One has nearly central collisions of two nearly equal ions. The experiment is 



carried out at various beam energies in the laboratory. We take the number of paticles to be 
85 (corresponding to central 40 Ar+ 45 Sc collisions). In a purely classical model the ground 
state has no kinetic energy at zero temperature so that the ground state energy per nucleon 
is —eN™' x / n where N™ x is the maximum number of nn bonds possible for particle number 
n = 85. Since N™ x is determined by geometry we can use experimental binding energy 
(rs 8.5 MeV/n) to fix the value of e. At temperature T the average energy per particle is 
l.hkT — eN nn /n where N nn , the average value of N nn is obtained from computer simulations. 
We can then write 

hT + e(NZ ax -Nnn)=e* (7.1) 
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For equal mass non-relativistic nuclear collisions we have e* = Eb eam /A where Eb eam is the 
beam energy per nucleon in the laboratory. There is an implicit assumption here that all 
available energy is converted to thermal energy. Thus the temperature is related to the 
beam energy. We can now fix the freeze-out density at different values, obtain an effective 
r at each beam energy and obtain points as in experiments |SIJ. This is shown in Fig. 10. 
However the fit with data is not good for any of the densities employed. 

This type of mapping between beam energy and temperature is not accurate. One of 
the sources of errors is the collective flow which is known to account for some fraction of 
the available energy. Better mapping could be expected where temperature is deduced from 



other experimental data p5| , pq| . In this approach, the tail of the proton cross-section is 
fitted by assuming that the proton has a Maxwell-Boltzmann distribution in a frame which 
is moving in the laboratory. We take this mapping from ||20|| . When the mapping from 
this phenomenological approach is used, the fit with the experimental data is quite good 
when the freeze-out density is taken to be 0.39p (see Fig. 11). What is also very pleasing 
is that the predictions for different freeze-out densities are also sufficiently different to be 
experimentally accessible. The combined study of S2 and r as a function of beam energy 
should be useful in determining the freeze-out density. 



VIII. HIGH TEMPERATURE CHARACTERISTICS 

Most of the attention in the present work has been focussed to temperatures that are 
close to what is believed to be the critical temperature of nuclear matter. Indeed, in the past 
the percolation model has mostly been used for mild to moderate excitation energies. At 
higher energies (i.e., Bevalac energies) other models have been used with moderate success. 
These models use approximations that are valid in the high temperature/low density limits. 
Chemical equilibrium between different species is assumed. We will call these models by a 
generic name, the thermodynamic model. In this model rather simple expressions for the 
average number of monomers (single nucleons= ni), dimers (clusters of two nucleons= 7T2), 



15 



three body clusters (723) etc. can be obtained. An early review can be found in J27 



We do not expect to find exact correspondence between these calculations and the present 
classical model at high temperature limits since we do not have quantum degeneracies. The 
clusters in our model consisting of attached cubical boxes in the lattice have degeneracies 
also. Even for moderate sized clusters these degenaracies require considerable effort to 
enumerate analytically. Nonetheless, many features seen in experiments are common in 
both the models. For example, in the thermodynamic model, n\, 7T2 and vary with 
temperature but in a way that T 3/2 n 2 /(ni) 2 remains constant. Here we have neglected the 
binding energy of the dimer with respect to kT. Constancy for this ratio is obtained in 
our present model also. Similarly nin 3 /(n 2 ) 2 is a constant in both the models. The so 
called coalescence relation ^p^-(2p) oc (^p^(p)) 2 is obeyed in both the models. Thus rather 
reasonable features emerge when the model is extrapolated to the high temperature side. 



IX. SUMMARY 

We have presented a model that has links with both mean field aspects and fragmentation 
of nuclei. We studied the equation of state under different approximations. The model 
shows a liquid-gas type phase transition. We also studied the nuclear fragmantaion and 
discussed the critical exponents near critical or percolation points. Some comparisons with 
experimental data were made. 

The purpose of this paper was to present the essentials of this model. The issues we 
addressed here are far from exclusive. Many other features may be further explored. The 
present model, we believe, is one step forward from the percolation model which was proved 
to be helpful for the analysis of experimental data. 



ACKNOWLEDGMENTS 

We would like to thank H. Guo, C. Gale, J. Kapusta and M. Grant for useful discussions. 
This work was supported in part by the Natural Sciences and Engineering Research Council 
of Canada and by the FCAR fund of the Quebec Government. 

16 



REFERENCES 

[1] M. W. Curtin, H. Toki and D. K. Scott, Phys. Lett. B123, 289 (1983). 

[2] H. Jaquaman, A. Z. Mekjian and L. Zamick, Phys. Rev. C27, 2782 (1983). 

[3] W. Bauer, D. R. Dean, U. Mosel and U. Post, Proc. seventh High Energy heavy ion 
study (GSI Darmstadt, 1984), GSI Report-85-10, p.701. 

[4] X. Campi and J. Debois, Proc. seventh High energy heavy ion study (GSI Darmstadt, 
1984), GSI Report-85-10, p.707. 

[5] D. Stauffer and A. Aharony, Introduction to percolation theory (Taylor and Francis, 
London, 1992). 

[6] W. Bauer, Phys. Rev. C38, 1297 (1988). 

[7] X. Campi, Phys. Lett. B208, 351 (1988). 

[8] K. Huang, Statistical Mechanics (John Wiley and Sons, Toronto, second edition, 1987). 
[9] D.H.E. Gross, L. Satpathy, T. Meng and M. Satpathy, Z. Physik C 309, 42 (1982). 
[10] P. J. Siemens, Nature 305, 410 (1983). 

[11] J. A. Lopez and P. J. Siemens, Nucl. Phys. A431, 728 (1984). 

[12] L.P.Csernai and J.I. Kapusta, Phys. Rept. 131, 223 (1986). 

[13] W. Bauer, G.F.Bertsch and S. Das Gupta, Phys. Rev. Lett. 58, 863 (1987). 

[14] J. Gallego, S. Das Gupta, C. Gale, S. J. Lee, C. Pruneau and S. Gilbert, Phys. Rev. 
C44, 463 (1991). 

[15] Ph. Chomaz, G. F. Burgio and J. Randrup, Phys. Lett. B254, 340 (1991). 

[16] D.H.E. Gross and B.H. Sa, Nucl. Phys. A437, 643 (1985); H.R. Jaqaman and D.H.E. 
Gross, Nucl. Phys. A524, 321 (1991). 

17 



[17] S. Ayik and C. Gregoire, Nucl. Phys. A513, 187 (1990). 

[18] K. C. Chase and A. Z. Mekjian, Phys. Rev. C49, 2164 (1994). 

[19] J. E. Finn, et al, Phys. Rev. Lett 49, 1321 (1982).. 

[20] T. Li, et al, Phys. Rev. C49, 1630 (1994). 

[21] T. Li, et al., Phys. Rev. Lett, bf 70, 1924 (1993) 

[22] A. Coniglio and W. Klein, J. Phys. A: 13, 2775 (1980). 

[23] M.E. Fisher, Physics, 3, 255 (1967). 

[24] R. B. Stinchcombe, in: Phase transitions and critical phenomena, vol. 7, eds. C. Domb 
and J. L. Lebowitz (academic press, London, 1983) 

[25] B.V. Jacak et.al, Phys. Rev. C35, 1751 (1987) 

[26] C.A. Ogilvie et. al., Phys. Rev. Lett. 67, 1214 (1991) 

[27] S. Das Gupta and A. Z. Mekjian, Phys. Rept. 72, 131 (1981) 



18 



FIGURE CAPTIONS 



Fig. 1 A square lattice is divided into blocks to illustrate the Bethe-Peierls approximations. 
See text for details. 

Fig. 2 The P — V diagrams in lattice gas model at temperatures T = 5, 10, 15 and 20 
(MeV). The solid curves are for the Bethe-Peierls approximation, and the dashed 
curves are for the Bragg- Williams approximation. In the calculation we used e = 9 
MeV. 

Fig. 3 The same as Fig. 2, but the dashed curves are for the mean field theory with Skyrme 
interactions. 

Fig. 4 The bond probability pb is plotted as a function of temperature. The solid curve is 
obtained from our model given in eq. (|4.2| ), and the dashed curve is the Coniglio and 



Klein model shown in eq. ( |4.3| ) 



Fig. 5 The mass yield distribution, Y(A) vs. A, for lattice N = 5 3 and n = 64, at tem- 
peratures T/Tc =0.5, 1.0, 1.5 and 2.0. Here Tq = 1.1275e is the thermal critical tem- 
perature. This value of Tq taken from EM is more accurate than the Bragg- Williams 



approximation. 

Fig. 6 The value of r is plotted as a function of temperature at different densities. 

Fig. 7 The value of S2 is plotted as a function of temperature at different densities. 

Fig. 8 The value of S2 is plotted as a function of multiplicity of fragments at different 
densities. 

Fig. 9 A max is plotted as a function of multiplicity of fragments at different densities. 

Fig. 10 The theoretical exponent r is compared with experimental data at different beam 
energies. The curves are obtained by using the temperature calculated from eqs. 
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J]). The solid circles are the corrected data taken from [21 1, the open circles are the 



uncorrected data taken from ||20|| , and the crosses are taken from |26| 



11 The same as Fig. 10, but the curves are obtained by using the temperature fitted 
from experimental data. 
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